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Abstract 

We study the charmless decays B — )• Khi] and B — )• Khi]' within the framework of QCD 
P^! factorization (QCDF) for K^ = K , K* , K^{Wi{)) and naive factorization for K^ = /s:|(1430). 

A . There are three distinct types of penguin contributions: (i) h — )• sqq — )• sr^q^ (ii) h — )• sss — )• sr]s, 

^ . and (iii) b — )• sqq — >• qKh, where r]q = {uu + dd)/\/2 and ry^ = ss. B ^ K^*'r]^ ' decays are 

dominated by type-II and type-Ill penguin contributions. The interference, constructive for Kr]' 
and K*rj and destructive for Krj and K*rj' , between type-II and type-Ill diagrams explains the 
pattern of T{B -^ Krj') > T{B -^ Krj) and T{B -^ K*r]') < T{B -^ K*7]). Within QCDF, the 
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VO . observed large rate of the Kr]' mode can be naturally explained without invoking flavor-singlet 

contributions or something exotic. The decay pattern for B — )• KQ{14:30)r]^ ' decays depends on 

ly-v , whether the scalar meson i^g (1430) is an excited state of k or a lowest-lying P-wave qq state. Hence, 

O , the experimental measurements of S — )• KQ{lA30)'r]^ ' can be used to explore the quark structure 

of iiro(1430). If iiro(1430) is a low-lying qq bound state, we find that KqTj has a rate slightly larger 
than KqT]' owing to the fact that the t]-?]' mixing angle in the rjq, rjs flavor basis is less than 45°, in 

k> \ agreement with experiment. Type-Ill penguin diagram does not contribute to i? — )• -R'l^^ ' under 

the factorization hypothesis and type-II diagram dominates. The ratio T{B — ^ K2f]')/T{B — t- -K^I^) 
is expected to be of order 2.5 as a consequence of (i) \ft\ > \f^\ and (ii) a destructive (constructive) 
interference between type-I and type-II penguin diagrams for K2I] (i^|ry'). However, the predicted 
rates of i? — t- i^|?7^ ' in naive factorization are too small by one order of magnitude and this issue 
remains to be resolved. There are two K^*'ri^ ' modes in which direct CP asymmetries have been 
measured with significance around Aa : Acp{K^r]) = — 0.37± 0.09 and Acp{K*^r]) = 0.19 it 0.05. 
In QCDF, power corrections from penguin annihilation which are needed to resolve CP puzzles 
in K~Tr^ and ir^ir^ modes will flip Acp{K^ri) into a wrong sign. We show that soft corrections 
to the color-suppressed tree amplitude 02 in conjunction with the the charm content of the rj will 
finally lead to Acp{K~r]) = — 0.15lg'28. Likewise, this power correction is needed to improve the 
prediction for Acp{K*^r]). 
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FIG. 1: Three different penguin contributions to B ^ KhV^'^ with Kh denoting K,K*,Kq (14:30) 
and 7^2(1430). Fig. 1(a) is induced by the penguin operators 03,5,7,9. 



TABLE I: Experimental branching fractions (in units of 10 ^) of i? — )• i^h^ with K^ = 
K,K*,K^{U30) and K^iUSO) taken from Q,Q. 

Kt] KrI K*r] K*r]' K^r] K^r]' K^r] K^r]' 

B+ 2.36 ±0.27 71.1 ±2.6 19.3 ±1.6 ^.'dt^l 15.8 ±3.1 5.2 ±2.1 9.1 ±3.0 28.0+||j 

5° l-12lo:28 66.1 ±3.1 15.9 ±1.0 3.8 ± 1.2 9.6 ±1.9 6.3 ± 1.6 9.6 ±2.1 13.7+|? 



I. INTRODUCTION 

Recently BaBar has measured charmless B decays with final states containing rj ov rj' [l|]. 
Comparing the first measurements of B ^ i^Q(1430)?7' and B — )• i^|(1430)T/' by BaBar with 
previous results of i? — > KQ{lA30)r] and B — > K2{1430)r] (see Table H]) clearly indicates that 
B{B -^ K^{U30)rj') < B{B -^ K^(1430)r/) and B{B -^ K^{U3Q)t]') > B{B -^ K^{U30)r]). It is 
well known that B{B -^ Krj') > B{B -^ Krj) and B{B -^ K*r]') < B{B -^ K*7]). The last two 
patterns can be understood as the interference between the dominant penguin amplitudes. 

For the rj and t]' particles, it is more convenient to consider the flavor states qq = {uu + dd)/^/2, 
ss and cc labeled by the rjq, rjg and 7?|?, respectively. Neglecting the small mixing with 7y|?, we write 

\ |r?) / \sm(p coscp J \ \r]s) ) 

where (\) = (39.3 ± 1.0)° [3|] is the rj — rj' mixing angle in the rjg and rjs flavor basis. Three different 
penguin contributions are depicted in Fig. [D (i) 6 — )• sqq — t- sr]q, (ii) b — )■ sss — >• sijs, and (iii) 
b — )• sqq — )• qKfi, corresponding to Figs. 1(a), 1(b) and 1(c), respectively. For B — )• K^*'r]^ > decays, 
the dominant penguin amplitudes arise from Figs. 1(b) and 1(c) governed by the parameters 
on{Kh.ris) and ai{r]qKh), respectively. Their expressions in terms of the effective Wilson coefficients 
04 and ae are summarized in Table [III 

It is clear that the interference between the B — t- Kr]q amplitude induced by the b — t- sqq penguin 
and the B — )■ Kr]s amplitude induced by 6 — )• sss is constructive for B -^ Krj' and destructive 
for B — )■ Krj. This explains the large rate of the former and the suppression of the latter [J]. 
For B — )• K*ri^ ' decays, it is the other way around. The sign difference between a4^{rjqK*) and 
a4^{K*r]s) explains why T(B — ;• K*r]) 3> T{B — )• K*r]'), recalling that 04 and qq are negative and 



TABLE II: The parameters a4,{Khris) and a^^irjqKh) with Kh = K,K*,K^{UZ{)) and ir|(1430). 

K K* K|^(1430) ^1(1430) 

ai{Khris) a4 + r5^=a6 a^-r'^aQ a^-r'^a^ a^-r'^^aQ 



the magnitude of the latter is larger than the former and that the chiral factor r^ to be defined 
below is of order unity for Hght mesons. 

The decay pattern for B — )• iTg (1430)7/*- ^ decays depends on whether irQ(1430) is an excited 
state of K, (or i^Q(800)) or a low lying P-wave qq state. Hence, the experimental measurements 
of -B — )• KQ{1430)ri^ ' can be used to explore the quark structure of the scalar meson Kq{1430). 
A detailed study in this work shows that T(B — )• K^r]) <^ T(B — )• KqTj') in the first scenario for 
E:^(1430) and T{B -^ K^r]) > T{B -^ K^t]') in the latter scenario. As for B -^ K^{U30)ri('^ 
decays, Fig. 1(c) does not make contribution owing to the vanishing decay constant of K2. Since 
the interference between Figs. 1(a) and 1(b) is constructive for K2ri' and destructive for -fr|?7 and 
since the decay constant /5 is larger than /* one will expect a larger rate for -fi'l^' than -ftTI'?- 

Recently we have studied the decays B — )• {K,K*)r]^ ' within the framework of QCD factoriza- 
tion (QCDF) [a, l6(. Here we shall present updated results with some discussions. Then in the rest 
of this work we will focus on i? — )• (/Cq(1430), iC|(1430))7/*^ ' decays and study their decay pattern. 

The layout of the present paper is as follows. In Sec. II we recapitulate the framework of 
QCD factorization. The we proceed to study B — )• {K,K*){rj,rj') decays in Sec. HI and B — >■ 
iiro(1430)(r/,r/') decays in Sec. IV. Since the QCDF approach for the i^|(1430)r/( > modes has not 
been developed, we reply on naive factorization to study the tensor meson production in Sec. V. 
Sec. VI comes to our conclusions. An appendix is devoted to the decay constants and matrix 
elements of the r/ and rj' mesons. 

II. QCD FACTORIZATION 

Within the framework of QCDF [8|, the effective Hamiltonian matrix elements are written in 
the form 

{M^M2\n,^\B) = ^ Y^ X^HM,M2\TA^+TBnB) , (2) 

V -^ p=u,c 

where Xp = VpbV*g with q = s,d, Ta describes contributions from naive factorization, vertex 
corrections, penguin contractions and spectator scattering expressed in terms of the flavor operators 
a^, while Tg contains annihilation topology amplitudes characterized by the annihilation operators 
b^. The explicit expressions of Ta and Tg can be found in [71, l8(. In practice, it is more convenient 
to express the decay amplitudes in terms of the flavor operators a^ and the annihilation operators 
Pf. Their relations to the coefficients of and 6f will be specified below. 

The expressions of -B — )• Khr]^"> decay amplitudes for K^ = K , K* , K^{143,Q) and i^|(1430) are 
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The order of the arguments of a^{MiM2) and /3f (M1M2), which are not shown exphcitly here, is 
consistent with the order of the arguments of the factorizable matrix elements x^^^i'*^^) given by 

X(BP„P,) ^ (P2|JlO)(Pi|j;|;B)=z/p,(m|-ml,jFo^^^(m|,J, 
XiBV,P) ^ {P\Ji^\0){V\Jl\B) =2fpmBPcAE''{ml), 

?BP/_,2\ 



Ml 



(4) 



X^BP^S) ^ (^S\J^\Q)(^p\j'\B) = fs{ml-ml)Ft [m-s,, 



{w?p) + m],b-{ml>)]e* p^p^B 



x(BP,T) _ (r|j^|o)(P|j;|s) = o, 

where fp,fv,fs are the decay constants of pseudoscalar, vector and scalar mesons, respectively, 
and k,b^,b- are B to tensor meson transition form factors defined in Eq. (j27p below. 

The flavor operators a^ and the annihilation operators /3f are related to the coefficients a^ and 
6fby 

ai(MiM2) = ai(MiM2), 

a2{MiM2) = a2{MiM2), 

j af (M1M2) - af (M1M2) for M1M2 = PP, VP, SP, TP 
\ af (M1M2) + ag(MiM2) for M1M2 = PV, PS, 
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al{MiM2) 
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a^(MiM2) + r^^2 ag(MiM2) for M1M2 = PP, PV, 
al{MiM2) - r^2 a|;(MiM2) for M1M2 = FP, 5P, PS, TP 

al{MiM2) - a?(MiM2) for M1M2 = PP, VP, SP, TP 

al{MiM2) + a?(MiM2) for M1M2 = PV, PS, 

a?o (Ml M2 ) + rf 2 aP (Ml M2 ) for Mi M2 = PP, PV, 

rf 2 al{MiM2) for M1M2 = VP , SP, PS, TP 



(5) 



<o(MiM2) 



and 



^f (M1M2) = U^IMiIMii>p 



(6) 



where the chiral factors r^^s are given by 

rPu\ = ?!!!£ ^ r'^s ^ ^P 

^ m;,(/i) fv ' ^ mbip,)ims-mg){n)' 

with the parameters /p and Zip being defined in the Appendix. 

The flavor operators af are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions [7|, |8| 

a^iM.M,) = (c. + ^) mM,) + ^ ^ [v^iM,) + ^H,iM,M,)] + PUM,), (8) 

where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), Cj are the Wilson co- 
efficients, Cp = {N^ — l)/{2Nc) with Nc = 3, M2 is the emitted meson and Mi shares the 
same spectator quark with the B meson. The quantities Vi{M2) account for vertex corrections, 
Hi{MiM2) for hard spectator interactions with a hard gluon exchange between the emitted meson 
and the spectator quark of the B meson and Pi{M2) for penguin contractions. The expression of 
the quantities Ni{M2) reads 

r 0, i = 6, 8 and M2 = y, 
Ni{M2) = \: ' (9) 

In Eq. ([3]), possible ffavor-singlet penguin annihilation contributions are denoted by /S^'s which 
will not be considered in this work. 

Power corrections in QCDF always involve troublesome endpoint divergences. For example, 
the annihilation amplitude has endpoint divergences even at twist-2 level and the hard spectator 
scattering diagram at twist-3 order is power suppressed and posses soft and collinear divergences 
arising from the soft spectator quark. Since the treatment of endpoint divergences is model depen- 
dent, subleading power corrections generally can be studied only in a phenomenological way. We 
shall follow [sj] to model the endpoint divergence X = /q dx/{l — x) in the annihilation and hard 
spectator scattering diagrams as 

Xa = In (^^ (1 + pAe''^^), Xh = In {^^ (1 + p^e^^^), (10) 

with Kh being a typical scale of order 500 MeV, and pa,h-, 4'A,h being the unknown real parameters. 
As pointed out in [6(], while the discrepancies between experiment and theory in the heavy quark 
limit for the rates of penguin-dominated two-body decays of B mesons and direct CP asymmetries of 
B^ — )• K^TT^, B^ —7- K^p^ and B^ — )• 7r"^7r~ are resolved by introducing power corrections coming 
from penguin annihilation, the signs of direct CP-violating effects in B~ — )• K~Tr^,B~^ — )• K~7] 
and B^ — )• vr'^vr'^ are ffipped to the wrong ones when confronted with experiment. These new B-CP 
puzzles in QCDF can be resolved by the subleading power corrections to the color-suppressed tree 
amplitudes due to spectator interactions and/or final-state interactions [9(] that not only reproduce 
correct signs for aforementioned CP asymmetries but also accommodate the observed B^ — )• vr'^Tr'' 
and p^TT^ rates simultaneously. Following [61], power corrections to the color-suppressed topology 
are parametrized as 

a2^a2il + pce''^n, (H) 



with the unknown parameters pc and 4>c to be inferred from experiment. 

For 5'-wave mesons, input parameters such as decay constants, form factors, quark masses, 
Wolfenstein parameters, hght-cone distribution amphtudes, power correction parameters can be 
found in [5|]. Input parameters for parity-even mesons such as iirQ(1430) and ir|(1430) will be 
specified later. For the renormalization scale of the decay amplitude, we choose /i = mb{mf,) = 4.2 
GeV. 

III. B -^ {K, K*){ri, 7]') DECAYS 
A. Branching fractions 

Details of the calculations in the framework of QCDF for all B — t- PP, VP decays can be 
found in [5|, la]. The updated results for the branching fractions and direct CP asymmetries in 
B —7- K^*>r]^ > decays are exhibited in Table Hill after correcting some minor errors in the previous 
computer codes. 

Numerically, Beneke and Neubert already obtained B{B~ — ;> K~if) ~ 0(50 x 10^^) in QCDF 
using the default values PA = Ph = m- Here we found similar results 57 x 10^® (53 x 10~^) 
with (without) the contributions from the "charm content" of the rj' . In the presence of penguin 
annihilation, we obtain B{B^ — )- K^rj') ~ 75 x 10^^ (67 x 10^^) with (without) the "charm 
content" contributions. Therefore, the observed large B — t- Krj' rates are naturally explained 
in QCDF without invoking, for example, significant flavor-singlet contributions or an enhanced 
hadronic matrix element (0|s75s|r/'). Data on i? — t- Kt] modes are also well accounted for by 
QCDF. 

The values of the parameters a4(i^'*''^q,s) and a4{r]qK^*') are given by 

a4{Kr]s) = -0.098 - 0.013i, a4{VqK) = -0.094 - 0.013z, 

a4iK*r]s) = 0.035 - 0.0009i, 04(7/^^*) = -0.036 - 0.0085i (12) 

The magnitude of ai^rjqK*) is smaller than ai{riqK) owing to the smallness of r^ ~ 0.36 compared 
to r^ ~ 1.45 at the scale p = 4.2 GeV, while the smallness of a4^{K*ris) relative to ai{Kris) is due 
to the destructive interference in the former. The sign difference between a/^{'qqK*) and ai{K*r]s) 
explains why T{B — )• K*r]) ^ T{B — )■ K*rj'). Although the rates of K*r]' and Kr] are comparable, 
B{B — )• K*vi) is much smaller than B{B — )■ Krj'). 

The QCDF prediction for the branching fraction of i? — ?• K*!]', of order 1.5 x 10^®, ^ is smaller 
than the predictions of pQCD and soft-collinear effective theory (SCET), but it is consistent with 
experiment within errors. The experimental values quoted in Table Hill are the BaBar measurements 
[l|. Belle obtained only the upper bounds: B{B~ -^ K*'r]') < 2.9 x 10"^ and B{B° -^ K*^r]') < 
2.6 X 10~^ uM- Therefore, although our central values are smaller than BaBar, they are consistent 
with Belle. It is very important to measure them to discriminate between various model predictions. 



1 The predictions BiB" -^ K*--q') = 2.2 x 10"^ and 6(5" -^ K*°7j') = 1.9 x lO'^ obtained by Beneke and 
Neubert [7| in the so-called "S4" scenario in which power corrections to penguin annihilation are taken 
into account are consistent with ours. 



TABLE III: Branching fractions (top; in units of 10~^) and direct CP asymmetries (bottom; in 
units of %) of -B — )• {K,K*)(ri,ri') decays obtained in various approaches. The pQCD results 
are taken from [lOt] for B — )• Kr]^ > with partial NLO corrections and from llj] for B -^ K*t]^ '. 
There are two solution sets with SCET predictions for decays involving r] and/or rj' [l2l. \ldt\. The 
theoretical errors correspond to the uncertainties due to the variation of (i) Gegenbauer moments, 
decay constants, quark masses, form factors, the A^ parameter for the B meson wave function, 
and (ii) pA,H, 4'A,H, respectively. 



Mode 



QCDF (this work) 



pQCD 



SCET 



Expt. [1, 2] 



B- 


-^ K rj 


9 9+1.7+1.1 
^•^-1.0-0.9 


q 9+1.2+2.7+1.1 
"J-^-0. 9-1. 2-1.0 


2.7 ±4.8 ±0.4 ±0.3 
2.3 ±4.5 ±0.4 ±0.3 


2.36 ±0.27 


B- 


^ K~i 


74 r:+57.9+25.6 
'^■^-25.3-19.0 


r-i n+13-5+11. 2+4.2 
O-'--'-'- 8.2- 6.2-3.5 


69.5 ±27.0 ±4.4 ±7.7 
69.3 ±26.0 ±7.1 ±6.3 


71.1 ±2.6 


B^ 


^ K°r] 


1 c + 1.4+0.9 
■■-•^-0.8-0.7 


9 1+0.8+2.3+1.0 
^•-■■-0.6-1.0-0.9 


2.4 ±4.4 ±0.2 ±0.3 
2.3 ±4.4 ±0.2 ±0.5 


1 19+0.30 
■■■•■■■^-0.28 


B^ 


^K^ri' 


70 q+54.1+24.2 


crn q+11.8+11. 1+4.5 
ou.o_ 8 2_ 6.2-2.7 


63.2 ±24.7 ±4.2 ±8.1 
62.2 ±23.7 ±5.5 ±7.2 


66.1 ±3.1 


B- 


^ K*-1] 


-, c- 0+8.2+9.6 
id.O_42_7.3 


09 1 q+0.26 
ZZ.iO_Q 27 


17 Q+5. 5+3.5 
■■■'•^-5.4-2.9 
IS R+4.5+2.5 


19.3 ±1.6 


B 


-^ K*-rl 


1 (.+2.1+3.7 
■■■•"-0.3-1.6 


6.38 ± 0.26 


A t:+6.6+0.9 
■^•O-3.9-0.8 

A + 5.3+0.8 

4.»_3.7_o.6 


A.8t\i " 


B^ 


^ K*^7] 


1 r 7+7.7+9.6 
-•■O-' -4.0-7.3 


QO Ql +0.28 
ZZ.01_o.29 


Ifi c;+4-l+2.3 
■■■"•0-4.3-2. 


15.9 ±1.0 


B^ 


-^ K*^7]' 


1 c+l. 8+3.5 
■■■•^-0.3-1.6 


q qcr+0.29 


A 1+6.2+0.9 
^•■■■-3.6-0.7 
A Q+4.7+0.7 
^•^-3.4-0.6 


q 1+1.0 b 
•^•■■■-0.9 


B- 


-^ K^r] 


-lA t;+10.3+15.5 
-'■^•'-'-26.0-10.7 


1 1 7+6.8+3.9+2.9 
■■^■■■•'-9.6-4.2-5.6 


33 ± 30 ± 7 ± 3 
-33 ±39 ±10 ±4 


-37 ±9 


B' 


^ K-T]' 


r, ^t;+0.69+1.20 
^•"^0-0.55-0.98 


f^ 9+1.2+1.3+1.3 
"•^-1.1-1.0-1.0 


-10±6±7±5 
0.7 ±0.5 ±0.2 ±0.9 


1.31}:? 


B^ 


^K^r, 


90 (<+ 9.8+12.6 
^<J-U_26.2-12.5 


19 7+4.1+3.2+3.2 
■■■^•'-4.1-1.5-6.7 


21±20±4±3 

-18 ±22 ±6 ±4 




50 


^K^rj' 


q n+0.6+0.8 
-^•^-0.5-0.8 


9 q+0.5+0.3+0.2 
z.O_o.4_o.6-0.1 


11±6±12±2 

-27±7±8±5 


5±5 


B- 


-^ K*-7] 


in 1+3.9+6.5 
-lU.l_37_7g 


-24.57t°:^? 


9 e+5.4+0.3 
^•"-5.5-0.3 
1 Q+3.4+0.1 

"■■■•^-3.6-0.1 


2±6 


B- 


-^ K*-7]' 


fjQ 7+ 6.5+27.9 
oy-' -38.6-49.5 


4mt\il 


9 7+27.4+0.4 
^•'-19.5-0.3 
9 R+26. 7+0.2 
^•"-32.9-0.2 


-26 ± 27 












B^ 


^K*^7] 


q 4+0.4+2.7 
'^•'i-0.4-2.4 


0.57 ±0.011 


1 1+2.3+0.1 
■"■•■■■-2.4-0.1 
n 7+1.2+0.1 
"•'-1.3-0.0 


19 ±5 


B^ 


-^ K*%' 


0+ 8.8+30.8 
°-°-10. 7-24.1 


-1.30 ±0.08 


q (-+ 8.9+1.3 
^•"-11.0-1.2 


2 ±23 



Q Q+6.2+0.9 
^•^-4.3-0.9 



''This is from the BaBar data 



"This is from the BaBar data 1|. Belle obtained an upper limit 2.9 x 10 ^ 1A\ 



Belle obtained aii upper limit 2.6 x 10 
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B. CP asymmetries 

There are two modes in which direct CP asymmetries have been measured with significance 
around Aa : AcpiK-rj) = -0.37 ± 0.09 and Acp{K*°r]) = 0.19 ± 0.05. It is crucial to understand 
them. Since the two penguin diagrams Figs. 1(b) and 1(c) contribute destructively to 5 — )• Krj 
due to the opposite sign of /^ and /^, the penguin amplitude is comparable in magnitude to the 
tree amplitude induced from b — )• usu, contrary to the decay B — )■ Kij' which is dominated by large 
penguin amplitudes. Consequently, a sizable direct CP asymmetry is expected in B~ — )• K~r] but 
not in K~7]' [15]. 

In the absence of any power corrections, it appears that the QCDF prediction Acp{K^r]) = 
— 0.233]'^o!i93 obtained in the leading \/mb expansion already agrees well with the data. ^ How- 
ever, this agreement is just an accident. Recall that when power corrections are turned off, the 
predicted CP asymmetries for the penguin-dominated modes K~-k^ , K*~7r~^, K~p~^, K~p^, and 
tree-dominated modes vr+vr", p^vr^ and p^vr"'" are wrong in signs when confronted with experiment 
[5J, Isl]. That is why it is important to consider the effects of power corrections step by step. The 
QCDF results in the heavy quark limit should not be considered as the final QCDF predictions to 
be compared with experiment. It turns out that the aforementioned wrong signs can be flipped 
into the right direction by the power corrections from penguin annihilation. 

However, a scrutiny of the QCDF predictions reveals more puzzles with respect to direct 
CP violation. While the signs of CP asymmetries in K~7r~^,K~ p^ modes etc., are flipped to 
the right ones in the presence of power corrections from penguin annihilation, the signs of Acp in 
B — )• K^TT^, K'r], TT~rj and B^ — t- vr^vr^, K*^r] will also get reversed in such a way that they 
disagree with experiment [5|,l6[. Specifically, Acp{K~7]) is found to be of order 0.11 in the presence 
of penguin annihilation and hence it has a wrong sign. These CP puzzles can be resolved by having 
soft corrections to the color-suppressed tree coefficient 02 so that 02 is large and complex. When pc 
and (f>c are turned on (see Eq. ([TT]) with pc = 1.3 and (ftc = —70° [5]), Acp{K~ri) will be reduced 
to —0.02 if there is no intrinsic charm content of the rj. Although the decay constant /^ ~ —2 MeV 
is much smaller than /^'« [see Eq. (1X12]) ]. its effect is CKM enhanced by VcbV*s/{VuhV*s). There- 
fore, the charm content of the r/ may have a significant impact on CP violation. Indeed, when /^ is 
turned on, Acp{K~rj) finally reaches the level of —15% with a sign in agreement with experiment. 
Hence, CP asymmetry in B~ — )■ K'-q is the place where the charm content of the ry plays a role. ^ 
The pQCD prediction is similar to QCDF. For comments on the pQCD calculation of Acp{K~rj), 
the reader is referred to |5|]. Note that while both QCDF and pQCD lead to a correct sign for 
Acp{K^ri), the predicted magnitude still falls short of the measurement —0.37 it 0.09. 

As for CP asymmetry in B^ — )• K*^7], we have Acp{K*^7]) = (3.ll^'g)% in the heavy quark limit. 
It is modified to (0.12l'l]^gg)% by penguin annihilation. Soft corrections to the color-suppressed 



^ Beneke and Neubert [Tl obtained Acp{K rj) = -O.lSOlggoo i^i their default predictions. 

■^ One of us (C.K.C.) has studied residual final-state interaction effects in charmless B decays |£(|. In 
this approach, Ac'p{K^ri) of order —0.27 can be induced through the decay B^ — > K^rj' followed by 
K^rj' — > K^rj rescattering via penguin diagrams. This rescattering mimics the effect of the charm content 
in the 77. 



{•iAtf 



tree amplitude is needed to improve the prediction and finally we obtain Acp{K 'q) — \'J-'^_2A) 
(Table lln|) . Unlike the K^r] mode, the charm content of the r/ here does not play an essential role 
for CP violation in K*^7]. We see from Table Hill that QCDF is in better agreement with experiment 
than pQCD and SCET, though it is still smaller than the data. 



IV. B -^ K^{U30)r]("^ DECAYS IN QCDF 



The hadronic charmless B decays into a scalar meson and a pseudoscalar meson have been stud- 
ied in QCDF in [16[. For scalar mesons above 1 GeV we have explored two possible scenarios in the 
QCD sum rule method, depending on whether the light scalars k, ao(980) and /o(980) are treated 
as the lowest lying qq states or four-quark particles: (i) In scenario 1, we treat k, ao(980), /o(980) as 
the lowest lying states, and iCg (1430), oq (1450), /o (1500) as the corresponding first excited states, 
respectively, and (ii) we assume in scenario 2 that i^g (1430), ao (1450), /o( 1500) are the lowest ly- 
ing resonances and the corresponding first excited states lie between (2.0 ~ 2.3) GeV. Scenario 
2 corresponds to the case that light scalar mesons are four-quark bound states, while all scalar 
mesons are made of two quarks in scenario 1. We found that scenario 2 is preferable. Indeed, 
attice calculations have confirmed that ao(1450) and Kq{1A30) are lowest-lying P-wave qq mesons 
18| , and suggested that a and k are 5- wave tetraquark mesonia [l8|, \l^ . 

Decay constants of scalar mesons are defined as 

(5'(p)l927/.ft|0) = fsPfi, {S\q2qi\0) = msfs- (13) 

the vector decay constant fs and the scale-dependent scalar decay constant fs are related by 
equations of motion 

ms 



fJ-sfs = fs, 



with ^5 



(14) 



m2(/i) -mi(//)' 

where 1712 and mi are the running current quark masses and ms is the scalar meson mass. In 
general, the twist-2 light-cone distribution amplitude (LCDA) of the scalar meson ^s has the form 



<^s{x,n) = fs6x{l 



l+f,sY.B^{^i)C'J\2x-l) 



m=l 



(15) 



where B^ are Gegenbauer moments and Ci 
we use 



.3/2 



are the Gegenbauer polynomials. For twist-3 LCDAs, 



<I>|(x) = /s, a>§(x) = /s6x(l 



(16) 



Since ^s = l/^o ^ 1 ^-nd even Gegenbauer coefficients are suppressed, it is clear that the LCDA of 
the scalar meson is dominated by the odd Gegenabuer moments. In contrast, the odd Gegenbauer 
moments vanish for the n and p mesons. The Gegenbauer moments Bi, B^ and the scalar decay 
constant fx* in scenarios 1 and 2 obtained using the QCD sum rule method 16|] are listed in Table 

lYl 



Form factors for B ^ S transitions are defined by 



{S{p')\A^\B{p)) 



m 



B 



mi. 



?BS, 



q, Fr{q') + 



m 



B 



TTIq 



BSi 



q,F,"'{q^) 



(17) 



TABLE IV: Gegenbauer moments Bi, B-^ and the scalar decay constant fx* (in units of MeV) in 
scenario 1 and scenario 2 at the scales /^ = 1 GeV and 2.1 GeV (shown in parentheses) obtained 
using the QCD sum rule method [16!]. 





Bi Bs 


Ik* 


scenario 1 
scenario 2 


0.58 ± 0.07 (0.39 ± 0.05) -1.20 ± 0.08 (-0.70 ± 0.05) -300 ± 30 (-370 ± 35) 
-0.57 ± 0.13 (-0.39 ± 0.09) -0.42 ± 0.22 (-0.25 ± 0.13) 445 ± 50 (550 ± 60) 


TABLE V: Values of a2,3,4 at /u = 4.2 GeV in scenarios 1 and 2. 




a2{KoVq,s) aaiK^riy^s) "4(^0%) "K^o^g) 


"4(^.^0) 


scenario 1 
scenario 2 


-0.27-0.08i 0.050 + 0.006i 0.025 + O.OOli 0.0001 - 0.0017i 
0.51-0.08i -0.047 + 0.006i 0.043 + 0.009i 0.018 - 0.002i 


0.50 + 0.09i 
0.45 + 0.02i 



where P^ = {p + p')^, q^i = (p — p')fi- As shown in [17|, a factor of {—i) is needed in i? — )• S* 
transition in order to obtain positive -B — )■ 5 form factors. This also can be checked from heavy 
quark symmetry [17(] . As a consequence, the factorizable amplitudes X^ ' ' and X^ ' ' defined 
in Eq. (HD are of opposite sign. In this work we shall use the form factors ^ 



-.BK, 



ni' 



0.21 



in scenario 1 and 



.BK* 



°iq' 



1 - 0.59(g2/m|) + Om{q^/m%) 



0.26 
1 - 0.44(g2/m|) + Om{q'^/m%) 



(18) 



(19) 



in scenario 2. They are obtained using the covariant light-front quark model 17l |. 
For the convenience of ensuing discussions, we would write 



A 



B^K^V 



A 



B^K'r,' 



^{BK*,r,q)(j^ j^ j^{BK*,-ns)(j^ j^ X^^'li'^fDC-^, 



(20) 



where Ci, C2 and C3 terms correspond to Figs. 1(a), 1(b) and 1(c), respectively. The expressions 
of Cj's in terms of the parameters a^ can be found in Eq. ([3]). As noticed in 16(, the pattern of Oj 
or Oj in i? — )• SP decays (see Table IV]) can be quite different from that in PP modes. For example, 
we have 



0.21 - 0.08z, 
-0.073 - O.Olli, 



a^{Kri, 



q,s) 



-0.010 - 0.0058i, 



(21) 



Oi2{Kriq^s 

aliKvq) = -0.073 - O.Olli, al{Kr]s) 

for Kr]^ > modes. Since the chiral factor ?'x " ~ ^^-"^ ^^ ^ 

one order of magnitude owing to the large mass of ErQ(1430), it follows that a^r/gi^Q) is much 



-0.098 - 0.013«, alir]gK) = -0.094 - 0.013i 



4.2 GeV is larger than r^ = 1.5 by 



16| , we have emphasized that the decay constant and the form factor for the excited 



In footnote 9 of 

state are of opposite sign. Hence, we do not agree with [l0|,|2l| on the sign of Fj^^" in scenario 1 



^BR- 



IO 



greater than aKi^gr/^) and a1{KQr]q). Likewise, the real part of a2,3(-f^o^(j,s) is greater than that 
of a2,3{K7]g^s) because the hard spectator term H{MiM2) 

with ^M and ^m being the twist-2 and twist-3 LCDAs of the meson M, respectively, is greatly 

enhanced for Mi = i('g(1430) due to the large chiral factor r^^'' . 

(') 
Because of the small vector decay constant of i^Q(1430), X^^'^i '^0) is suppressed relative to 

(') * (') 

X^^^h^i ) and X^^^o''^!! ). However, C3 gains a large enhancement from aKr/g/fg). As a result, 

the amplitude of Fig. 1(c) is comparable to that of Fig. 1(a). The decay pattern of i? — t- i^g'? 

depends on whether the scalar meson i^g(1430) is an excited state of k or a lowest-lying P-wave 

qq state. In scenario 1, we have (in units of GeV'^) 

XiBK*'^,) = -0.60, x(^^o-'?=) = 0.61, ^(^'''•■^o) = -0.15, 

XiBK*^i^) = -0.49, x(^^o.'7D = _o.75, x(^<'^o) = -0.12. (23) 

Consequently, Fig. 1(b) interferes destructively (constructively) with Figs. 1(a) and 1(c) for KqT] 
[KqI]'). This leads to B{B -> KqT]') > B{B -> KqI]) (see Table EI]), which is in sharp disagreement 
with experiment. 

The decay pattern in scenario 2 is quite different. Because of the large magnitude of a3{KQr]g^s) 
and the large cancellation between as{KQr]s) and a4{KQrjs) in C2, B — )• KqTj^ ' decays are dominated 
by the contributions from Figs. 1(a) and 1(c), contrary to the B — )■ K^*>r]^ > decays which are 
governed by Figs. 1(b) and 1(c). Numerically, we obtain 

Ci = -0.098 + O.Olli, x(^^o.'?9) = _o.74, x^^'^'^'-^o) = 0.22, 
C3 = 0.455 + O.OlTf, X^^^^o'V'g) = _o.60, X^-^^'-^o) = 0.18 . (24) 

Therefore, the penguin diagrams Figs. 1(a) and 1(c) contribute constructively to both K^rj and 
K^T]' with comparable magnitudes. Since X^^^^'V,) /x^^K^^v'^i) = x^^'ii'^S) /X^^'^'-^'^o) = cot cp^ 
1.23, it is clear that AB-^K*q/AB~^K*ri' ^ cot0 and hence B — ;■ KqT] should have a rate larger than 
B — )• KqT]' in scenario 2 as the mixing angle (p is less than 45°. 

The QCDF results for branching fractions and CP asymmetries of i? — )• -ftrg(1430)77^ ' obtained 
in two different scenarios are shown in Table IVII where the central values are for pA = pc = and 
0A = 0c = and the ranges < pa,h < 0.5 and < (J)a,h < 27r have been considered for the 
estimate of uncertainties. It is evident that scenario 2 is in better agreement with experiment than 
scenario 1. This implies that the scalar meson Kg (1430) is a lowest lying rather than excited qq 
state. 

Several remarks are in order, (i) If X^^-^o''??) and X^^^i'^o) are of the same sign, the interference 
between Ci and C3 terms will become destructive in scenario 2 and yield too small rates for both 
K^T] and K^rj' . (ii) A recent pQCD calculation indicates that B{B -^ K^r]')/B{B -^ K^rf) ^ 2 
in scenario 1 and B{B — )■ KqT]') ~ 75 x 10~^ in scenario 2. Neither of them is consistent with 



.BK* 



experiment, (iii) Different estimates of form factors at g = were obtained recently: Fg " (0) 



-0.34 



21[, -0.44 [la] in scenario 1 and Fg^^°(0) ~ 0.60 [21[, 0.76 ^ in scenario 2. If these form 



1 X 10 ^ in scenario 2 for Fg " (0) = 0.60 as an example. This implies that a small form factor 



factors are used in QCDF calculations, we will have B{B — )• KqT]) ~ 22 x 10 ^ » B{B — > KqT]' 
1 X 10"^ in scenario 2 for Fq^^{0) = 0.6 
for B — >• i^g(1430) transition is preferable. 
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TABLE VI: Branching fractions (top; in units of 10 ^) and direct CP asymmetries (bottom; in 
units of %) oi B ^ K^{U?>Q){r],r]') decays obtained in QCDF (this work) and pQCD [20] in 
scenario 1 (first entry) and scenario 2 (second entry). 



Mode 



QCDF (this work) 



pQCD 



Expt. [1, 2j 



S-^ A'*"(1430)?7 
B- ^ K*- {WiQ)r]' 



5° -^ i^o*°(1430)T/ 



50 -^ K^^{Um)r]' 



r, c:+1.7+11.9 
'-'•^-0.4- 0.5 
-10 q+6.2+48.6 
^"J-'-'-a. 6-10.1 



11. 8^ 



-5.3+0.3+1.1+2.5 
-3.5-0.4-1.2-2.3 
33 g+13. 5+1. 1+7.7+8.2 



27.8: 

9.6:1 



1-20.8+14.1 

-10.2-17.6 

12.9+10.0 

6.1- 9.4 



9.0-1.1-7.0-7.3 
+1.6+3.1+4.0+4.5 
-0.5-2.8-3.6-4.1 
rjj 5+15.8+6.2+21.0+18.0 



21.6 



r, 4+1.5+11.8 

^•^-0.3- 0.4 
19 «+5. 7+50.1 
-■-^•0-3.4- 9.3 



9.1 



28.4l^^; 



26.6 

8.7 



19.7+17.3 
9.8-16.5 
+11.7+7.5 
- 5.5-3.1 



22.0 



10.8-5.8-16.5-16.1 

-4.4+0.0+1.1+2.0 

-2.8-0.1-1.1-1.8 

6+1.4+6.4+6.9 

8-1.4-5.9- 
+1.6+3.2+3.94 



6.2 
4.6 
-0.5-3.6-3.0-4.2 
Y/l 9+15.0+6.4+20.5+17.2 
'^•^-10.3-5.7-16.2-15.5 



15.8 ±3.1 

5.2 ±2.1 
9.6 ±1.9 

6.3 ±1.6 



B- ^ K*Q- {Um)ri 



B- ^ A'*-(1430)r?' 



^0 ^ i^*o^^43Q)^ 



^0 ^ i^*0(i43Q)^/ 



8.3 



-21.1+11.7 



- 5.4 

1.7 



-0.3 

1.2 

14.5 

2.2 



11.7 
0.2+2.4 
-0.2-2.4 
+0.2+0.4 
0.2-0.4 
+0.7+1.7 
-0.6-1.7 

36.3+20.4 
10.8-20.4 
+0.2+3.1 
0.2-3.1 



_n n+0.1+0.3 
^•^-0.1-0.3 
-, -1+0.5+1.6 
^•-^-0.5-1.6 



5 ±13 

6 ±20 
6 ±13 

-19 ±17 



V. B ^ K2*(1430)7?(') IN NAIVE FACTORIZATION 

The observed J^ = 2+ tensor mesons /2(1270), /:^(1525), 02 (1320) and ATI (1430) form an 
SU(3) 1 ^^2 nonet. The qq content for isodoublet and isovector tensor resonances is obvious. The 
polarization tensor e^j^ of a ^P2 tensor meson with J^^ = 2'^~^ satisfies the relations 



0, 



^^u 



^^lU 



where p^ is the momentum of the tensor meson. Therefore, 

(0|(y - A)^\Tie,p)) = ae^^p" + be^p^ = 







(25) 



(26) 



and hence the decay constant of the tensor meson vanishes identically; that is, the tensor meson 
cannot be produced from the V — A current. This means that Fig. 1(c) does not contribute to 
B —7- K2ri^ ' under the factorization hypothesis. 



The general expression for the B 
{T{e,pT)\{V-A)^\B{pB)) = 



T transition has the form 



22| 



ih{q^)e^,vp.e*'"'pBo.{pB+PTY{pB - PtY + k{q^)e%PB 



+ h+{q^)elpp%p^,,{pB + Pt)^ + h^{q^)elpp%p^^{pB - Pt)^- (27) 



The form factors h, k, 6+ and 5_ have been calculated in the ISGW quark model [22| and its 
improved version, the ISGW2 model J23l |. They have also been computed in the covariant light- 
front (CLE) quark model 17| and are listed in Table IVTIl for form factors fitted to a 3-parameter 
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TABLE VII: Parameters in the form factors of i? — )• -fC|(1430) transitions in the parametrization 
of Eq. (|28p . as obtained by fitting to the covariant hght-front model IT]]. The form factor k is 
dimensionless, while k, b^ and 6_ are in units of GeV~^. The numbers in parentheses are the form 
factors at q^ = obtained using the ISGW2 model 23 1. 



F 


F(0) 


a 


b 


F 


F(0) 


a 


b 




0.008 (0.016) 
-0.006 (-0.006) 


2.17 
1.96 


2.22 
1.79 




0.015 (0.293) 
0.002 (0.0063) 


-3.70 
0.38 


1.78 
0.92 



form 



F(^i 



F(0) 



(28) 



1 — a{q'^ /m?^) + b{q'^ /m?^Y ' 

Notice that when q^ increases, h{(f') , 6+(g^) and b-{q^) increases more rapidly in the CLE model 
than in the ISGW2 model and that the form factor k in both models is quite different. 
The decay amplitude of i? — )• TP always has the generic expression 



A{B ^ TP) = e* p>^ M{B ^ TP) . 



The decay rate is given by 

T(B 



TP) 



Pc 



fmB\' 



\M{B -^TP)\'^ 



(29) 



(30) 



127rm|, \mT J 

where Pc is the magnitude of the 3-momentum of either final-state meson in the rest frame of the 
B meson. Since the QCDE approach for B — t- TP has not been developed, we will reply on naive 
factorization to make estimates; that is, we will not include vertex, hard spectator and penguin 
corrections in Eq. ([8]) for the calculation of ai . 

The predictions of i? — t- i^|(1430)ry' •* in various models are shown in Table IVnll We see that 
the predicted rates in naive factorization are too small by one order of magnitude. ^ It was found 



m 



24| that B{B — )• K2'ri')/B{B — )• -^^2^) ~ 45, a prediction not borne out by experiment. This is 



ascribed to the fact that the matrix element 

(r?'|s75s|0) 



mt 



2m 



'f?i'^ 



(31) 



used in j24l ] is erroneous as it does not have the correct chiral behavior in the chiral limit; the 
SU(3)-singlet rji acquires a mass of the QCD anomaly which does not vanish in the chiral limit. 
Applying the anomalous equation of motion Eq. (|A7p and neglecting the masses of u and d quarks, 
one gets [27|, bd] 



(r?(')|s75s|0) 



m 



.(') 



2?7lo 



/,,(') 



1 



fQ 



(32) 



Although the form factor k{q ) is very different in the CLF and ISGW2 models, the magnitude oi F{q ) 



■2\ - 



k{q^) 



[m. 



ml)b+{q^) + qH.iq^ 



-0.15 in the former and 0.14 in the latter for q^ in the range 



between m^ and m"^, is about the same. Consequently, the predicted rates of B 
independent of the form-factor models. 



if|77^ ^ are basically 
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TABLE VIII: Branching fractions (in units of 10 ^) of i? — )• i^|(1430)(r/,r/') decays obtained in 
various approaches. The predictions of 2J] are for 1/N^^ = 0.3 . 



Mode 



This work 



KLO [24] 



MQ [25] 



SV [26] 



Expt. [1, 2j 



B- ^K;-{U30)t] 


1.1 


0.031 


1.19 


0.01 


9.1 ±3.0 


B- ^ K;- {U30)ri' 


2.7 


1.4 


2.70 


0.007 


28.0tli 


BO ^ Kf{U30)rj 


1.0 


0.029 


1.09 


0.01 


9.6 ±2.1 


gO -^ Kf(U30)r]' 


2.5 


1.3 


2.46 


0.006 


13.7^3-2 



Since /^,/v2 is about half of f^, , this means that the rj' matrix element of the pseudoscalar density 
is reduced roughly by a factor of 2. On the contrary, the t] matrix element is enhanced owing 
to the opposite sign of /^ and /^. A rigorous result without neglecting light quark masses is 
{"n I^TS'SJO) = —W f^,J{2ms) [cf. Eq. ()A4p ]. Since the magnitude of /i^ is numerically close to /iL, 
see Eq. (]A12p . it becomes clear that the rate of A'lr/ is comparable to but slightly smaller than 
the K2V1' one. Our results agree with 251] as we use the same form factors derived from the CLF 
model. ^ The predictions of 261] are too small as the authors only considered the tree diagram 
effects and neglect the important penguin contributions. 



A slightly large rate of A'lV over i^|?7 comes from two sources: (i) \ft\ > J/^J and (ii) a 
destructive (constructive) interference between Figs. 1(a) and 1(b) for -fCI'? (-^IV)- It is expected 
that T(B — )• K2'i]')/T{B — t- A'I?]) ~ 2.5 . Since the predicted absolute rates are too small by one 
order of magnitude, ^ the discrepancy between theory and experiment may call for the study of (i) 
NLO effects from vertex, hard spectator scattering and penguin corrections, (ii) power corrections 
from penguin annihilation, and (iii) improved estimate of i? — )■ T form factors. The detailed 
investigation will be left to a future work. 



VI. CONCLUSIONS 

We have studied the decays B — )■ Khr] and B — )• Khi]' within the framework of QCDF for 
Kfi = K,K*,Kq{143,Qi) and naive factorization for K^ = i^2(1430). There are three different types 
of penguin contributions: (i) b — )• sqq — )• srjq, (ii) h — )• sss — )• srjs, and (iii) h — )• sqq — )• qK^, 
corresponding to Figs. 1(a), 1(b) and 1(c), respectively. Our conclusions are as follows: 

• i? — )• K^*>r]^ > decays are dominated by type-II and type-Ill penguin contributions. The 
interference, constructive for Kr]' and K*'q and destructive for Kr] and K*r]', between Figs. 
1(b) and 1(c) explains the pattern of r{B -^ Krj') > r{B -^ Krf) and T{B -^ K*r]') <C 



4]) for the rj^ ' matrix elements of 



^ The only difference between our work and 25] is that we use Eq. 
pseudoscalar densities, while the authors of 25] use Eq. ([5^ . 

^ The same issue also occurs in the D -^ TP system l29|. The predicted branching fractions based on fac- 
torization are at least two orders of magnitude smaller than data, even for decays free of weak annihilation 
contributions. We cannot find possible sources of rate enhancement. 
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T(B —7- K*rj). Within QCDF, the observed large rate of the Kr]' mode can be naturally 
explained without invoking flavor-singlet contributions or something exotic. The predicted 
central values of the decay rates for K*^ri' and K*^r]' are smaller than BaBar but consistent 
with Belle. 

There are two K^*'ri^ ' modes in which direct CP asymmetries have been measured with 
significance around 4cr : Acp{K-r]) = -0.37 ± 0.09 and Acp{K*^r]) = 0.19 ± 0.05 . In 
QCDF, power corrections from penguin annihilation which are needed to resolve CP puzzles 
in K^tt'^ and vr'^vr" modes will flip Acp{K~rj) into a wrong sign. We show that soft 
corrections to the color-suppressed tree amplitude 02 in conjunction with the the charm 
content of the rj will finally lead to Acp{K~rj) = — 0.15lo!28- ^^^ corrections to 02 are 
also needed to improve the prediction for Acp{K*^r]). The QCDF prediction is in better 
agreement with experiment than pQCD and SCET. 

The decay pattern of i? — t- i^Q(1430)ry( •* depends on whether i^Q(1430) is an excited state 
of K or a lowest-lying P-wave qq state. If iirQ(1430) is an excited state of gg, one will have a 
destructive (constructive) interference of Fig. 1(b) with Figs. 1(a) and 1(c) for K^r] {Kq-tj'). 
This leads to B{B -^ K^r]') > B{B -^ K^r]). If K^(1430) is made of the lowest-lying qq, 
we found that Figs. 1(a) and 1(c) interfere constructively and that A{B — t- KQ'q)/A{B — )• 
KqT]') ~ cot(j) with (j) being the rj-rj' mixing angle in the rjq,ris flavor basis. Hence, KqT] has 
a rate slightly larger than K^rj' owing to the fact that is less than 45°. The agreement 
of the latter scenario with experiment indicates that the scalar meson i^g(1430) is indeed a 
bound state of the low-lying qq state in P-wave. 

Fig. 1(c) does not contribute to i? — )■ i^^^' •* under the factorization hypothesis and Fig. 
1(b) dominates. The ratio T{B — t- K2rj')/T{B — t- i^^^) is expected to be of order 2.5 as 
a consequence of (i) \f5\ > \f^\ and (ii) a destructive (constructive) interference between 
Figs. 1(a) and 1(b) for Xl^ (-R'lV)- However, the predicted rates of S — )■ -f^l^^ ' in naive 
factorization are too small by one order of magnitude and this issue remains to be resolved. 
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Appendix A: The rj — rj' system 

Decay constants /^(/j, /'^(/) and /^^/j are deflned by 

(01^7^7591/^) = i^/J,)P^, (0|s7^75s|/)) = i/^(/)P^, (0|c7^75c|?7(')) = i/^(/)P/., (Al) 
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while the widely studied decay constants fq and fs are defined as 

mri5q\r]g) = -^Up^ , {0\sYl5s\7is) = ifsP^ . (A2) 



The ansatz made by Feldmann, Kroll and Stech (FKS) |3| is that the decay constants in the quark 
flavor basis follow the same pattern oi r] — rj' mixing given in Eq. ([T]) 



Empirically, this ansatz works very well [3|]. Theoretically, it has been shown recently that this 
assumption can be justified in the large- A^c approach 31|. 
Consider the matrix elements of pseudoscalar densities [30| 

2m,(0|g75(z|r?(')) = ^/^Jo. 2m,(0|s75s|r?(')) = i/i^^,,, (A4) 

one can define the parameters hq and hs in analogue to fq and fs 

2mg{0\qj5q\rjg) = — hg , 2m^(0|s75s|?7s) = ihs , (A5) 

and relate them to /i^'L by the same FKS ansatz 



h?) ^ri \ ( cos <f> — sin (f) \ I hq 
K' K' / V ^"^ ^ cos(/) J \ hs 



(A6) 



Using the equations of motion 

47r 

one can express all non-perturbative parameters in terms of the decay constants /g,/s and the 
mixing angle 0: 



d^^is^l^s) = 2im, S75S + -^ G^, G^" , (A7) 



^g = /g(?7i^cos^ (/) + m^, sin^</)) — \/2/s(m^/ — m^)sin(/)cos(/), 

hs = fsi'm-n'^^^ (j) + m sm. (j)) ^(m / — tti ) sin(/)cos0. (A8) 

v2 

For numerical calculations we shall use those parameters determined from a fit to experimental 
data PI 

/g = (1.07 ±0.02)/^, /, = (1.34 ±0.06)/^, </. = 39.3° ± 1.0°. (A9) 

The masses of r/g and ry^ read (3| 

Jg Jg 47r /^ 47r 

m^^ = l(0|m,.-i75s|r?.) + -^(0|^GG|r/,)«2m|,-m2 + -i(0|^GG|r?,), (AlO) 

where contributions to their masses from the gluonic anomaly have been included. We shall use 
the parameters extracted from a phenomenological fit [3(]: 

^(0|^GG'|r?g)= 0.265 ±0.010, _J0M^^^ = k^ (All) 
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Numerically, 



/i« = 0.0013 GeV^ /i^ = -0.0555 GeV^ h'^, = 0.0011 GeV^ /i^, = 0.068 GeV^ 

/^ = 109 MeV, /^ = -111 MeV, /^, = 89 MeV, f^, = 136 MeV, 

/^ = -2.3MeV, /^/ = -5.8MeV, m^, = 741MeV, m^, =802MeV, (A12) 



where we have used the perturbative result [32] 
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